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ON  Q-MATRICES,  CENTROIDS  AND  SIMPLOTOPES 


by 

Richard  W.  Cottle 
and 

Rabe  von  Randow 


1.  Background .  Our  concern  in  this  paper  is  the  question:  What  is 
the  class,  Q  ,  of  all  real  n  x  n  matrices  M  such  that  for  every  n- 
vector  q  the  linear  complementarity  problem  (LCP) 

Iw  -  Hz  =  q 
v  ^  0  ,  z  >_  0 
(  w,  z>  =0 

has  a  solution? 

Several  researchers  have  already  attacked  the  problem  of  character¬ 
izing  such  matrices,  and  -  in  a  sense  -  solved  it.  See  [1],  [6],  [7],  [9]. 
Many  more  studies  have  identified  subclasses  of  Q  as  well  as  methods 
for  solving  the  linear  complementarity  problems  formed  thereby.  But  the 
evidence  to  date  suggests  that  a  "useful"  characterization  of  Q  will 
not  be  attained  by  accumulating  larger  and  larger  classes  of  matrices 
contained  in  it. 

Lately,  a  frontal  assault  on  the  geometry  of  the  problem  has  begun. 

In  the  vanguard  of  this  attack  are  L.  M.  Kelly  and  L.  T.  Watson  to  whose 
papers  [5],  [6]  this  research  owes  much.  In  particular,  we  adopt  their 
reformulation  of  the  problem  as  one  of  covering  the  sphere  by  spherical 
simplices.  We  now  endeavor  to  make  this  statement  more  precise. 

The  LCP  can  be  cast  in  what  looks  like  (but  ultimately  need  not  be) 
more  general  terms:  Given  real  n  *  n  matrices  A  and  B  and  a  real 
,  find  n-vectors  w  and  z  such  that 

Aw  +  Bz  =  c  (1  ') 

w  _>  0  ,  z>0  (2') 

<  w,  z)  =  0  O') 

♦ 


(1) 

(2) 

(3) 


n-vector  c 


•t 
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Let  us  designate  the  system  (1')  -  (3')  by  the  triple  (A,  B,  c)  .  This 
means  (1)  -  (3)  is  the  triple  (I,  -M,  q)  .  The  question  about  Q-matrices 
now  becomes:  What  are  necessary  and  sufficient  conditions  on  A  and  B 
to  ensure  that  (A,  B,  c)  has  at  least  one  solution  for  every  c  ? 

Since  (2')  and  (3')  imply 

wiZi  =  0  i  =  1 . n  ,  (4) 

attention  turns  to  what  are  loosely  called  complementary  submatrices  of 
[A,  B]  and  the  cones  they  span.  An  n  x  n  matrix  C  is  said  to  be  a 
complementary  submatrix  of  [A,  B]  if  and  only  if 

ci  €fV  V  '  (5) 

In  (5),  as  elsewhere,  At  ,  Bi  and  C  denote  the  ith  columns  of  A  , 

B  and  C  ,  respectively.  In  principle,  there  are  2°  complementary 
submatrices  of  [A,  B]  ;  but,  if  At  =  Bi  for  some  i  ,  there  will  be 
fewer  than  2  distinct  complementary  submatrices  of  [A,  B]  .  At  the 

moment,  such  duplications  are  of  no  consequence.  We  let 

comp  [A,  B]  (6) 

denote  all  n  x  n  matrices  C  formed  in  accordance  with  (5). 

Associated  with  each  C  €  comp  [A,  B]  is  the  so-called  complemen¬ 
tary  cone 


pos  C  -  {c:  c  =  Cx  ,  x  >  0}  . 


(7) 
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Solving  (A,  B,  c)  can  be  regarded  as  a  matter  of  finding  C  G  comp  [A,  B] 
such  that  c  G pos  C  .  In  fact,  the  union  of  the  complementary  cones 
(relative  to  [A,  B]  )  is  precisely  the  set  of  n-vectors  c  for  which 
(A,  B,  c)  has  a  solution.  What  interests  us  is  finding  necessary  and 
sufficient  conditions  on  A  and  B  for  which 

U  {pos  C:  C  G  comp  [A,  B]}  “  E  ,  (8) 

n 

where  En  denotes  Euclidean  n-space.  This  formulation  of  the  problem 
was  first  used  by  Samelson,  Thrall  and  Wesler  (10] and  later  by  Murty  [8]. 

A  nonsingular  matrix  C  G  comp  [A,  B]  is  called  a  complementary  basis 
in  [A,  B]  .  Notice  that  if  A  is  nonsingular,  the  problem  (A,  B,  c) 
can  be  transformed  to  (I,  A  ^B,  A  ^c)  which  more  closely  resembles 
(I,  -M,  q)  . 

Normally,  the  matrix  M  in  (I,  -M,  q)  is  called  nondegenerate  if 
and  only  if  every  principal  minor  of  M  is  nonzero.  It  is  easy  to  see 
that  M  is  nondegenerate  if  and  only  if  every  member  of  comp  [I,  -M] 
is  nonsingular.  With  this  in  mind,  we  call  [A,  B]  nondegenerate  if  and 
only  if 


det  C  ^  0  for  every  C  €  comp  [A,  B]  .  (9) 

Otherwise  (A,  B]  is  degenerate. 

From  the  formulation  of  (A,  B,  c)  it  is  evident  that  the  columns 
Ai’  Bi  (1  =  •••»  n)  and  c  can  be  positively  scaled  without  affecting 

the  solvability  of  the  problem.  Now,  since  (A,  B,  c)  is  trivial  when 
c  *  0  ,  we  may  assume 


u 


HI  =  i 


Furthermore,  it  can  be  shown  [ 12]  that  (8)  cannot  hold  if  any  column  of 
[A,  B]  is  0  .  Hence  in  the  present  investigation,  we  shall  assume 


I  ^-jJI  II  ®-jjl  i  =  1,  . . . ,  n  . 


This  means  that  c  and  all  the  columns  of  [A,  B]  can  be  regarded  as 
points  on  the  unit  (n-l)-sphere 


s"-1  -  (x  6En:  |„|  -  1) 


with  center  at  the  origin,  0  . 

For  the  present,  we  strengthen  the  assumption  (11)  by  assuming  that 

[A,  B]  is  nondegenerate.  Although  this  is  a  strong  assumption,  it  is 

far  from  being  powerful  enough  to  imply  much  vis-A-vis  the  existence  of 

solutions.  Under  the  present  nondegeneracy  assumption,  if 

C  =  [C  ,  ...,  C  ]  £  comp  [A,  B]  ,  then  the  vectors  C, ,  ....  C  are 
-in  In 

linearly  independent  and  (by  virtue  of  their  uniform  length)  can  be  con¬ 
sidered  as  the  vertices  of  a  spherical  (n-l)-simplex,  o(C)  .  Note  that 


o (c)  =  pos  c  ns 


In  line  with  (8),  the  question  can  be  put  as  follows:  What  are  neces¬ 
sary  and  sufficient  conditions  on  A  and  B  to  ensure  that 


U  fo(C):  C  e  comp  [A,  B] }  »  S  ? 
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We  want  to  know  when  the  collection  of  spherical  (n-l)-simplices  o( C) 
covers  the  (n-l)-sphere  S°  ^  . 

Following  Kelly  and  Watson  [6]  to  whom  this  interpretation  of  the 
problem  is  due,  we  shall  call  comp  [A,  B]  a  Q-arrangement  if  and  only 
if  (13)  holds;  and,  having  taken  the  liberty  of  using  the  term  "non- 
degenerate"  in  reference  to  [A,  B]  ,  we  go  a  step  further  and  call 
[A,  B]  a  Q-matrix  if  comp  [A,  B]  is  a  Q-arrangement.  It  is  clear 
that  with  A  being  nonsingular,  [A,  B]  is  a  Q-matrix  in  the  new  sense 
if  and  only  if  -A  '''B  is  a  Q-matrix  in  the  original  sense. 

We  propose  to  investigate  the  characterization  of  Q-matrices  accord¬ 
ing  to  the  formulation  given  by  (13).  Indeed,  we  state  a  geometrical 
conjecture  on  necessary  and  sufficient  conditions  for  (13)  to  hold.  We 
show  that  the  condition  is  only  necessary  in  general  and  that  it  is  also 
sufficient  in  the  case  n  =  2  . 

In  order  to  state  the  conjecture,  we  continue  with  the  development 
of  the  geometrical  structures  involved.  Let 

<*  £v  I’h11 

be  the  unit  ball  in  E  with  center  at  0  .  Note  that  Sn  ^  is  the 
-  n 

boundary  of  .'-S?n  .  Now  for  each  C  £  comp  [A,  B]  define 

t(C)  =  pos  C  n  -^n  .  (14) 


We  call  such  sets  spherical  sectors.  Our  nondegeneracy  assumption  implies 
that  the  t (C)  are  n-dimensional,  and  of  course  they  are  compact  and 
convex.  Thus  they  are  convex  bodies  [2j.  Moreover,  their  union  is  ,jBn 
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if  and  only  if  [A,  B]  is  a  Q-matrix. 

Now  we  regard  ,£n  as  a  convex  body  with  homogeneous  distribution 
of  mass.  For  each  C  £  comp  [A,  B]  let 


x (C)  =  centroid  (center  of  mass)  of  t(C)  . 


(15) 


If  v(C)  denotes  the  volume  of  x(C),  then 


x(C) 


1 

v(C) 


/ 

T(C) 


xdx  . 


(16) 


For  the  moment,  we  say  no  more  about  how  these  centroids  are  actually 
found.  Finally,  we  define 


X[A,  B]  =  conv  (x(C) :  C  £  comp  [A,  B] }  .  (17) 

That  is,  X[A,  B]  is  the  convex  hull  of  the  centroids  of  the  spherical 
sectors  relative  to  [A,  B]  .  At  last  we  can  state  the  conjecture  we 
want  to  study. 

Conjecture.  Let  [A,  B]  be  nondegenerate.  Then 

[A,  B]  is  a  Q-matrix  if  and  only  if  0  £  int  X[A,  B]  .  (18) 

As  mentioned  earlier,  we  shall  prove  that  (18)  holds  when  n  =2  and 
also  prove  that  for  all  n  >_  2  the  necessity  (0  £  int  X[A,  B])  must 
hold.  We  shall  also  investigate  the  converse  and  prove  it  is  false  for 
n  >  3  . 


T 
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2.  Proof  of  the  conjecture  for  the  case  n  =  2  .  It  must  be  conceded 
that  a  criterion  for  Q-matrices  is  not  very  interesting  in  the  case  where 
n  =  2  .  There  are  several  reasons  for  this.  Among  them  are  the  following 
arguments : 

1°  One  can  draw  the  complementary  cones  pos  C  and  see  whether  their 
union  is  E7  . 

2°  There  already  exist  finite  numerical  tests  for  whether  a  2*2 

matrix  M  belongs  to  Q  .  (See  [1],  [6].)  The  transformation  in¬ 
dicated  in  Section  1  coupled  with  such  a  test  would  handle  the  "more 
general"  case  of  [A,  B]  . 

To  these  may  be  added  another  objection: 

3°  The  result  -  being  a  theorem  in  plane  geometry  -  is  probably  already 
known.  (The  authors  have  not  yet  found  a  reference,  however.) 

These  points  notwithstanding,  there  is  a  simple  rationale  for  es¬ 
tablishing  (18)  for  n  =  2  :  It  makes  the  ideas  more  concrete. 

In  the  2-dimensional  case,  the  "spherical"  1-simplices  o(C)  are 
just  closed  arcs  on  the  unit  circle.  The  nondegeneracy  assumption  means 
that  for  each  C  =  [C  ,  C7]  £  comp  [A,  B]  ,  and  C?  are  neither 

identical  nor  antipodal  points.  The  length  of  n(C)  must  be  positive 
and  less  than  tt  .  A  "spherical"  sector  r(C)  is  just  an  ordinary 
circular  sector,  and  its  centroid  lies  on  the  line  segment  between  the 
center  of  the  circle  and  the  midpoint  of  o(C)  . 

Theorem  1.  Tf  A  and  B  are  real  2x2  matrices  such  that  [A,  B] 
is  nondegenerate  then 

[A,  B]  is  a  O-matrix  if  and  only  if  0  G  int  X[A,  B]  . 

T 

-  ... - - - -  - iW  tfl - MBWI I  '  ' . — 


Proof .  Suppose  [A,  B]  is  a  Q-matrix.  If  0  9  int  X[A,  B]  ,  then 
X[A,  B]  must  lie  within  a  half-disk,  the  region  bounded  by  a  semi-circle 
and  the  corresponding  diameter,  ST  .  Let  N  denote  the  midpoint  of 
the  other  semi-circle.  Thus  ON  j  ST  and  ST  strictly  separates  N 
from  X[A,  B]  .  Now  since  [/.,  B]  is  assumed  to  be  a  Q-matrix,  there 
exists  a  matrix  C  G  comp  [A,  B]  such  that  N  Ea(C)  •  Let  R  be 
the  midpoint  of  o(C)  .  Then  x(C)  ,  the  centroid  of  >  (C)  ,  lies  on 
the  radius  OR  ,  and  as  x(C)  belongs  to  X[A,  B]  ,  ST  strictly  separ¬ 
ates  x(C)  and  thus  R  from  N  .  But  this  implies  the  length  of  o(C) 
is  greater  than  ti  which  is  impossible.  Thus  0  E  int  X[A,  B]  .  (See 
Figure  1.) 

Conversely,  suppose  0  G  int  X[A,  B]  and  [A,  B]  is  not  a  Q- 
matrix.  Then  there  exists  a  point  N  on  S^  which  does  not  belong  to 
any  of  the  arcs  a(C)  .  Since  0  E  int  X[A,  B]  ,  there  cannot  be  an  en¬ 
tire  semi-circle  of  such  "uncovered"  points.  Hence  any  uncovered  point 
such  as  N  must  belong  to  an  open  minor  arc  whose  endpoints  are  A^ 
and  B^  ;  say  i  =  1  .  The  antipodal  points  corresponding  to  A^  and 

are  -A^  and  -B^  ,  respectively.  Since  N  is  uncovered,  it  follows 
that  both  A^  and  B^  must  belong  to  the  open  arc  from  -A^  to  -B^  . 
But  this  forces  all  the  centroids  x(C)  to  lie  in  the  same  half-disk, 
namely  the  one  determined  by  the  bisector  ST  of  the  central  angle 
-A^OB^  and  containing  -A^  and  -B^  .  (See  Figure  2.)  This  is  impos¬ 
sible  if  0  G  int  X[A,  B]  .  D 

Th  ■'rem  1  establishes  the  conjecture  (18)  for  the  limited  case  of 
n  =  2  .  We  mention  in  passing  that  one  can  readily  see  how  much  the 
nondegeneracy  assumption  can  be  relaxed.  When  n  =  2  one  needs  at  least 
three  distinct  nonsingular  members  of  comp  [A,  B]  . 


r  - 
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3.  Proof  of  necessity  in  the  general  case.  As  we  shall  demonstrate 
in  this  section,  the  separation  argument  used  in  the  first  half  of  the 
proof  of  Theorem  1  lends  itself  to  higher-dimensional  generalization. 

For  any  n  _>  2  ,  suppose  [A,  B]  is  a  nondegenerate*  Q-matrix  and 
0  £  int  X[A,  B]  .  Then  by  standard  separation  arguments  there  exists 
a  hyperplane  H  which  separates  0  from  X[A,  B]  ,  i.e.,  0  and 
X[A,  B]  do  not  lie  in  the  same  open  half-space  with  respect  to  H  . 

In  fact,  we  may  choose  H  to  satisfy 

H  =  fx:  (  d,  x)  =0)  (19) 

fidfl  =  1  (20) 

<  d ,  y>  0  for  all  y  £  X[A,  B]  .  (21) 


This  means  H  strictly  separates  d  from  X[A,  B]  .  Now  since  [A,  B] 
s  supposed  to  be  a  Q-matrix,  d  £  o( C)  for  some  C  £  comp  [A,  B]  . 

Let  t(C)  be  the  corresponding  spherical  sector,  and  let  x(C)  be  its 
centroid.  Since  x(C)  £  X[A,  B]  the  inequality 


<d,  x(C)>  >  0  (22) 

would  clearly  contradict  (21).  As  we  shall  show,  this  is  precisely  what 


*The  full  strength  of  the  nondegeneracy  assumption  is  not  needed  here. 

What  we  want  to  do  is  restrict  attention  to  the  nonsingular  C  £  comp  (A,  B]  . 
If  [A,  B)  is  a  0-matrix,  then  it  follows  that  Sn-1  is  covered  by  the 
spherical  (n-l)-simpl ices  a(C)  corresponding  to  these  nonsingular 
C  £  comp  [A,  B)  .  In  this  more  natural  setting  X[A,  B]  would  be  the 
convex  hull  of  the  corresponding  x(C)  . 


a. 


I, nl  iiHiir»)«fc 


is  implied  by  the  hypotheses  stated  above.  To  this  end,  we  first  prove 


Lemma  1.  The  inequality  (22)  holds  for  all  d  £o(C)  if  and  only  if 

<C.,  x(C)>  >0  i  =  1,  n  .  (23) 

Proof.  The  necessity  is  obvious  since  the  C.  are  points  of  o(C)  . 

For  the  sufficiency,  note  that  for  any  d  6o(C)  we  have 

n  n 

d  =  }’  X.C.  ,  X.  >  0  ,  y  A.  >  0  .  (24) 

,'.ii  i  —  .  .  l 

i=l  i=l 

Combining  (23)  and  (24)  we  obtain  (22).  □ 

Thus,  our  task  comes  down  to  establishing  (23).  In  geometrical 
terms,  (23)  says  that  the  centroid  of  a  spherical  sector  makes  an  acute 
angle  with  the  vertices  C^,  ...,  .  (Recall  that  in  the  case  of 

n  =  2  ,  this  is  obvious  since  x(C)  lies  on  the  bisector  of  a  central 
angle  of  less  than  rr  radians.) 

At  this  point,  our  work  requires  only  the  nonsingularity  of 
C  =  [C^,  C  ]  ,  the  associated  sets  cr(C)  ,  t(C)  ,  and  of  course 

x(C)  .  The  inequality  (23)  is  a  geometrical  statement  independent  of  the 
complementarity  setting.  It  should  also  be  pointed  out  that  the  fact  we 
are  about  to  establish  may  already  be  known,  but  the  authors  have  not 
succeeded  in  locating  it  in  the  literature. 

The  proof  will  be  made  somewhat  tidier  if  we  settle  a  few  minor 


; 
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points  beforehand.  First,  it  is  clear  that  a  proof  of  (23)  for  i  =  n 
which  makes  no  special  assumptions  about  can  be  adapted  to  prove 

x(C))  >  0  for  all  i  .  Second,  the  proof  makes  use  of  an  appro¬ 
priately  defined  wedge  containing  t(C)  .  The  wedge  is  split  into  parts, 
each  of  which  is  a  convex  body  in  its  own  right.  In  this  way,  we  relate 
the  centroids  of  the  parts  to  the  centroid  of  the  entire  wedge.  This 
device  makes  use  of 

Lemma  2.  Let  K  be  a  homogeneous  convex  body  with  volume  v(K)  and 

centroid  x(K)  .  If  K  is  partitioned  into  m  convex  bodies 

K, ,  ...,  K  with  volumes  v(K,),  ...,  v(K  )  and  centroids 
I  m  I  m 

x(K  ),  ...,  x (K  )  respectively,  then 
l  m 

_  m  v(K  )  _ 

x(K)  =  £  xO^)  .  (25) 

Proof.  This  follows  from  (16).  □ 

m 

Note  that  since  v(K)  =  £  v(K.)  ,  the  equation  (25)  means  that 

i=l  1  _ 

x(K)  is  a  convex  combination  of  the  x(K^)  ,  i  =  1,  ...»  m  . 

We  are  now  ready  to  state  and  prove  the  result. 

Theorem  2.  (c  ,  x(C))  >  0  . 


Proof .  We  begin  by  defining  the  (n-l)-dimensional  subspace 
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H  =  lx  £E  :  (C  ,  x)  =  01  . 

0  n  n 

The  hyperplane  Is  just  the  orthogonal  complement  of  the  1-dimensional 

space  generated  by  .  It  determines  two  half-spaces 

H^=  (x  £E  :  (c  ,  x)  >0) 

0  n  n  — 

and 

H~  =  {x  e  E  :  <  C  ,  x>  <  0}  . 

U  n  n  — 

There  are  three  cases  determined  by  the  way  in  which  the  vertices 
C^,  are  situated  with  respect  to  these  half-spaces.  To 

describe  these  cases  efficiently,  we  introduce  another  harmless  assump¬ 
tion.  Indeed,  we  may  assume  without  loss  of  generality  that  the  ver¬ 
tices  are  labeled  in  such  a  way  that 

<C±,  Cn>  £<Ci+1,  Cn>  i  =  1,  ....  n-2  .  (26) 

Case  1.  (C^,  C  ^  •  This  case  is  trivial  since  x(C)  is  a  positive 

linear  combination  of  the  (1  £  i  £  n)  .  In  view  of  the  fact  that 
<Cn,  Cn>  =  1  ,  we  have  the  required  inequality  ( C  ,  x(C)).  >  0  . 

Case  2.  < C, ,  C  >  <0  and  (C  ..  C  )  <  0  .  Here  we  define  a  set  of 
1  n  n-1  n  — 

n  more  distinct  (n-l)-dimensional  subspaces  (hyperplanes).  For 


i  *  1 . n  let 


H  •  {x  G  E  :  x 
1  n 


II 

l  Vi  • 

j=i  j  j  1 


Each  hyperplane  can  be  put  In  the  form: 


Hi  “  {x  G  En:  <Pi’  **  =  0}  1  “  1 . . 


where  the  vector  p^  satisfies 


<  Pj.  C±)  >  0  i  =  1 . n  . 


Again  each  of  these  hyperplanes  gives  rise  to  a  pair  of  half-spaces: 


H±  =  (x  €  En:  <  pit  x>  _>  0}  i  =  1,  ....  n 


Hi  °  {x  E  V  <Pi*  -  0}  i  -  1,  ....  n  . 


Note  that  the  (i  =  1,  n)  are  the  bounding  hyperplanes  of 

pos  C  .  Moreover, 


t:  =  t(C)  =  H*  . 

i=l  1 


Using  the  half-spaces  H^,  we  define  the  spherical  wedge 
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n-1 

w  =  .*?"  n  iT  . 

i=l 


(We  note  in  passing  that  t  =  Wn  H  .)  The  wedge  W  is  n-dimensional . 

n 

Let  denote  its  centroid.  It  is  easy  to  verify  the  following 

important  facts  about  W  .  First,  W  contains  the  diameter  of  jf/1 

passing  through  C  (and  the  antipodal  point  -C  ).  Second,  the  hyper- 
n  n 

plane  splits  W  into  two  "congruent"  parts.  Third  (and  really 

crucial) , 


<Cn’  V 


The  latter  follows  from  Lemma  2  and  the  remark  about  congruence. 
Now  notice  that  if  we  let 


n-l 

y:  n  hT  ri  h"  =  w  n  h* 

i-1  1  n 


then  we  can  write 


W  =  t  U  y  . 


The  set  y  is  also  a  convex  body.  Let  x  =  x(C)  ,  and  let  x  denote 

t  y 

the  centroid  of  y  .  Then 


v(t)  —  .  v(y)  — 

X  3  — - —  -  v  4-  — v 

w  v(W)  T  v(W)  Y 


(29) 
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where  v(t)  ,  v(y)  ,  and  v(W)  are  the  volumes  of  t  ,  y  ,  and  W  , 
respectively.  Since  y  is  a  subset  of  Hq  (by  the  assumptions  of  the 
present  case)  we  have 

<cn»  *  0  •  (30) 

The  desired  result  now  follows  from  (27),  (29),  and  (30).  (Figure  3 
depicts  the  situation  for  the  case  where  n  =  3  .) 

Case  3.  (C^,  cn^  <  0  <  (C^^,  cn>  •  This  hypothesis  means  that  the 

open  half-spaces  relative  to  HQ  each  contain  at  least  one  point  from 
the  set  {C^,  ...»  C^_^}  .  Under  these  circumstances,  the  idea  is  to 
reduce  the  proof  to  the  two  preceding  cases  by  introducing  some  new  point 
and  corresponding  convex  bodies.  (See  Figure  4.)  A  spherical  sector 
to  which  Case  1  applies  will  be  called  a  set  of  type  1,  and  a  sector 
to  which  Case  2  applies  will  be  called  a  set  of  type  2.* 

Of  particular  interest  is 

{Pij:  Pij  =  H0  n°(Ci’  cj>  •  ci  e  int  H~  ,  C  €  int  ,  j  /  n}  (31) 


*In  Figure  4,  the  sector  with  vertices  0,  P.2»  C_,  C.  is  of  type  1  where 
as  that  with  vertices  0,  C^,  P^2»  C^  is  of  type  2. 
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where  0(0^,  )  denotes  the  (1-dimensional)  edge  of  a(C)  with  end¬ 
points  and  C.  .  Each  P  defined  in  (31)  is  the  central  projec¬ 
tion  onto  Sn  1  of  a  point  C^.  lying  in  Hq  n on  the  line  segment 

between  C,  and  C.  .  The  points  C  ,  . ..,  C  ,  are  vertices  of  an 
i  J  i  n-1 

ordinary  (n-2)-simplex  A^,  . ..,  ^)  lying  in  .  What  we  want 

to  do  is  use  the  points  C..  ,  . .  . ,  C  ,  and  the  P..  to  define  a  parti- 

1  n-1  ij 

tioning  of  o( C.  ,  ....  C  ,)  into  (n-2)-simplices  such  that  all  the 
1  n-i 

vertices  of  each  individual  (n- 2) -simp lex  lie  either  in  Hq  or  Hq  . 

It  will  suffice  to  do  this  to  A(C  ,  C  )  using  only  C.,  ,  .  . .  ,  C  , 

i  n-1  1  n-1 

and  the  C_  as  vertices  of  the  (n-2)-simplices  in  the  partitioning. 

We  claim  that  a  partitioning  of  this  sort  exists.  (See  the  Appen¬ 
dix.)  By  adjoining  0  and  to  each  of  the  (n-2)-simplices  (of  the 

partition)  that  lie  in  H+  ,  we  obtain  spherical  sectors  of  type  1.  The 
centroid  of  each  such  sector  belongs  to  int  Hq  .  Similarly,  by 
adjoining  0  and  to  each  of  the  other  (n-2)-simplices  (of  the 

partition)  we  obtain  spherical  sectors  of  type  2.  (Each  of  these 
is  the  portion  of  a  wedge  lying  in  H  q  . )  By  Case  2,  the  centroid 
of  such  a  set  belongs  to  int  Hq  .  Assembling  all  the  aforementioned 
spherical  sectors  (on  both  sides  of  Hq  )  we  obtain  the  given  sector 
t  .  Now,  by  Lemma  2, 


x(C)  =  xt  €  int  Hq 


and  this  completes  the  proof.  □ 
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4.  Disproof  of  the  sufficiency  In  the  general  case.  In  the  preceding 
section,  we  remarked  that  the  nondegeneracy  assumption  could  be  dropped 
provided  that  X[A,  B]  were  properly  redefined.  Let  us  assume  this  is 
done. 

It  is  clear  that  the  property  0  G  int  X[A,  B]  is  not  adversely 
affected  by  slight  perturbation  of  some  columns  of  [A,  B]  .  This 
means  that  if  there  exists  a  Q-matrix,  [A,  B]  ,  having  an  arbitrarily 
small  perturbation  [A,  B]  which  is  not  a  Q-matrix,  then  the  property 
0  €  int  X[A,  B]  cannot  be  sufficient  for  [A,  B]  to  be  a  Q-matrix. 

The  existence  of  just  such  a  matrix  has  been  discovered  by  Watson 
[13].  Let  [A,  B]  =  [I,  -M]  where 

1  -1  4 

M  =  4  -3  1 

10  0 

The  matrix  M  is  degenerate  since  two  of  its  principal  minors  are  0  . 
To  see  that  M  is  a  Q-matrix  (equivalently  that  comp  [I,  -M]  is 
a  Q-arrangement)  takes  a  bit  of  work  which  we  shall  omit.  The  important 
point  though  is  that  for  e  >  0  ,  the  matrix 

1  -1  4 

M(e)  =4-31 
1  e  0 

is  not  in  Q  since  (I,  -M(e),  q(e))  has  no  solution  for 

0 

q(t)  =  1  , 

-e/4 


r  >  0 
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This  example  shows  that  for  n  _>  3  ,  0  £  int  X[A,  B]  does  not 

imply  that  [A,  B]  is  a  Q-matrix.  In  a  sense  this  is  easy  to  "explain." 
In  this  example,  there  are  six  centroids  x(C)  where  C  €  comp  [I,  -M] 
and  det  C  ^  0  .  The  origin  must  belong  to  the  convex  hull  of  four 
centroids  (though  not  necessarily  to  its  interior).  At  any  rate,  it 
seems  as  though  there  could  be  enough  degrees  of  freedom  to  insure  that 
0  €  int  X[A,  B]  and  still  arrange  the  vectors  A^  ,  B^  ,  i  =  1,  2,  3 
in  such  a  way  that  [A,  B]  is  not  a  Q-matrix. 


4UHIMWIIU.U..I  UCIWUWWVi 
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Thus 


C,  G  int  H  and  C  ,  €  int  . 
1  (J  n-X  0 


We  may  partition  the  vertices  C, ,  . . .  ,  C  ,  as  follows: 

1  n-1 


V"  =  { C± :  <Cn,  Ct>  <  0}  =  {C^, 


C,  ) 


,0 


-  { C± :  <Cn,  Ct>  0}  -  {Ck+1,  ....  C5 


V+  ’  {V  <C„’  Cl>  >  01  '  (Ct+l . Cn-1» 


with  the  stipulations  1  _<  k  <  S.+1  <_  n-1  and  the  understanding  that  V 

may  be  empty  in  which  case  k  =  2.  . 

Setting  C  aside  temporarily,  we  concentrate  on  the  face 
n 

A(Cp  Cn_^)  which  is  "split"  into  two  pieces  by  .  That  is, 

A(C^,  ....  C  has  vertices  lying  in  each  of  the  open  half-spaces 

determined  by  HQ  .  Therefore,  if  (C^,  C  ^  )  f=V  *  V+ ,  the  edge  ACC^,  C^) 
meets  in  a  point  which  we  denote  by  C  .  We  now  define  the  set 


G:  =  {CA  eEn:  C±.  =  HQ  Ga^,  C^)  ,  (C^  C  )  €  v"  x  V+}  , 


The  cardinality  of  G  is  k(n-l-2)  . 
Our  aim  is  to  partition 


A  (C^,  •••»  *  =  H0  ^  MC^,  •  ••» 


and 


4‘<ci . C„-1):  ■  »0  n4<ci . C„-l> 


J 
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into  (n-2)-simplices  whose  vertices  all  belong  to  V+  U  v*  U  G  and 
V  U  v'  U  g  ,  respectively,  with  pairwise  intersections  of  lower 
dimensions. 

It  is  easy  to  see  that  A  (C^,  C^_^)  is  (n-2) -dimensional , 

for  it  is  convex  and  contains  the  (n-2)-simplex 

A(CI’  (V  Ck+1’  Cit’  Cl,l+1*  Cl,n-1)  ‘  The  extreme 

points  of  A  (C, ,  ....  C  ,)  are  precisely  the  set  V  U  U  G  . 

I  n-1 

(Note  that  the  elements  of  V  U  ,  being  extreme  points  from  the 

start,  are  still  extreme.  The  points  of  G  must  also  be  extreme  since 

they  lie  on  distinct  edges  of  A(C^,  ...,  C^_^)  ,  each  having  one  of  its 

endpoints  not  in  A  (C,  ,  ...,  C  .)  .) 

i  n-l 

Our  proof  of  the  existence  of  the  partitioning  described  above  is 
based  on  the  special  case  where  V*"*  =  0  .  But  this  assumption  is  not 
so  restrictive,  for  it  can  be  achieved  by  a  small  perturbation  of 
(which  tilts  Hq  slightly).  Once  the  partition  of  the  perturbed  set 
is  obtained,  reversing  the  process  causes  some  of  the  simplices  in  the 
partitioning  to  disappear;  what  remains  is  still  a  partitioning  of  the 
required  type. 

Assume  =  <!>  .  Then  A  (C,  ,  ....  C  , )  is  a  frustum  of  a  simplex. 

1  n-l  - 

Indeed  the  hyperplane  H  separates  the  vertices  C. ,  . . . ,  C  into 

u  1  n-l 

two  sets:  V  (of  cardinality  k  )  and  V+  (of  cardinality  n-l-k  ). 

Following  Sommerville  (11,  p.  1031,  we  refer  to  A  (C, ,  C  ,)  as 

i  n-l 

a  frustum  of  type  (k| n-l-k)  ,  the  other  frustum 

4+<Cl . -  »o  na(cr  ••••  C„-l> 

being  of  type  (n-l-k|k)  .  The  first  entry  in  the  symbol  (*|*)  rep¬ 
resents  the  number  of  vertices  of  the  original  simplex  in  the  frustum  we 
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"keep"  and  the  second  entry  Is  the  number  of  vertices  "cut  off"  by  Hq  . 
Sommerville  shows  that  these  two  frusta  are  isomorphic  to  simplotopes, 

that  is,  Cartesian  products  of  simplices.  In  particular,  a  frustum  of 

type  (r|s)  is  isomorphic  to  A  *  A  .  To  illustrate,  consider 

r~ i  s 

the  frustum  of  type  ( 3 | 1 )  which  arises  when  a  plane  section  cuts  one 
vertex  off  of  a  tetrahedron.  The  frustum  is  clearly  isomorphic  to  a 
prism,  i.e.,  A^  x  A^  . 

The  objects  we  actually  partition  will  be  simplotopes,  Af  x  Ag  . 
We  may  assume  that  r  _<  s  since  the  order  of  the  factors  is  not  essen¬ 
tial.  However,  we  do  want  to  be  sure  that  once  this  convention  is  ex- 

tablished,  then  A  *  Ag  and  A^  x  Ag,  are  isomorphic  if  and  only  if 

r  =  r'  and  s  =  s'  .  Now  if  A  x  A  and  A  ,  x  A  ,  are  isomorphic, 

r  s  r  s 

they  have  the  same  number  of  vertices  and  faces  (of  maximal  dimension). 
The  number  of  vertices  is 

(r  +  l)(s  +  1)  -  <r'  +  D(s'  +  1)  • 


Thus 


rs+r+s  =  r's'+r'  +  s'  . 


The  number  of  faces  is 


r+s+2  =  r'+s'  +  2  . 


Therefore,  Ar  x  Ag  is  isomorphic  to  Af  .  x  Ag  ,  if  and  only  if 
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rs  =  r 's ' 
r+s=r'  +  s' 
r  j_  s 
r  '  <_  s  ' 

These  conditions  hold  if  and  only  if  r  =  r'  and  s  =  s"  •  See  Figure  7. 

■  r+s  '  *  r't'S  * 

Theorem  3.  The  simplotope  A  *  A  can  be  partitioned  into  |  r  I  =  t  s  i 
-  r  s  1 

simplices  of  dimension  r+s  • 


Proof.  The  r  +  s  +  2  faces  of  A  *  A  are  of  two  types.  Those  of 
-  r  s 

the  first  type  (of  which  there  are  r  +  1  )  are  of  the  form  A  ,  *  A 

r  -1  s 

Those  of  the  second  type  (of  which  there  are  s  +  1  )  are  of  the  form 


A^  x  A^  ^  .  The  generation  of  A^  x  Ag  can  be  performed  as  follows. 

(See  Sommerville  [11,  p.  113].)  Let  Pn  denote  a  vertex  of  A  .  Each 

U  r 

(r-l)-face  of  A  generates  an  (r+s-l)-face  of  A  x  a  when  P„  is 
r  r  s  U 

moved  all  over  Aq  .  This  gives  rise  to  r  +  1  (r+s-l)-faces  of  the 


form  A  ,  x  a  .  Also,  A  generates  an  (r+3-l)-face  of  A  x  a  when 
r-1  s  r  r  s 

P  is  moved  all  over  a  (s-l)-face  of  A  .  This  gives  rise  to  s  +  1 

s 

(r+s-l)-f aces  of  the  form  Af  x  A  ^  .  It  is  clear  from  the  way  the 


product  is  generated  that  the  intersection  of  an  (r-ks-l)-face  of  the  first 


type  and  an  (r+s-l)-face  of  the  second  type  is  an  (r+s-2)-face  of  type 


A  .  x  a  .  .  Furthermore,  the  intersection  of  two  (r-l)-faces  of  A 
i  —  i  s—  x  r 

is  a  (r-2)-face  of  the  form  2  •  Hence  two  (r+s-l)-faces  of  the  first 
(second)  type  have  as  intersection  an  (r+s-2)-face  of  the  form  A^  2  x 
(Ar  x  As_2)  . 

Now  let  F^  be  an  (r+s-l)-face  of  the  first  type  and  let  F^  be  a 

(r+s-l)-face  of  the  second  type.  Then  there  exists  exactly  one  vertex 

U  of  Af  x  which  does  not  lie  on  F^  U  F2  •  By  the  above,  F^  H  F2 

is  a  (r+s-2)-face  F,  of  type  A  ,  x  a  ,  .  Now 

3  r-1  s-1 
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(  r(s+l)  J 

has  /  (r+l)s  \  vertices  . 

(  rS  1 

Hence  F^  U  F2  has  r(s+l)  +  (r+l)s  -  rs  =  (r+l)(s+l)  -  1  vertices, 
exactly  one  less  than  Ar  x  As  *  Using  U  as  apex,  form  the  pyramids 
Fj  and  F 2  on  F  and  F 2  ,  respectively.  These  have  a  common 
(r+s-l)-face  F^  ,  the  pyramid  over  F^:  =  F^  O f2  .  Now  the  hyper¬ 
plane  generated  by  F,  cuts  A  *  A  into  two  pieces  K.  and  K»  . 

■j  x  s  1  2 

It  follows  readily  from  the  definitions  of  F^  ,  F2  ’  F3  an^  t^ie 

convexity  of  these  pyramids  and  of  K,  ,  K.  and  A  *  A  that 

1  l  r  s 


K1  = 


K2  " 


and  K, 


nv 


Hence 


F,  U  F„ 


A  x  A  . 

r  s 


The  partitioning  process  is  now  clear.  Given  Af  *  Ag  there  are 
faces  F^  ,  F 2  of  the  two  types  and  a  unique  vertex  V  contained  in 
neither  one  and  corresponding  pyramids  F^  ,  F^  such  that 
Fl  U  F2  =  &r  x  .  Taking  partitions  of  F^  and  F2  in  the  obvious 
inductive  manner  induces  simplicial  partitionings  of  F^  and  F2  . 
Assembling  chese  we  get  the  required  partitioning  of  A^  *  A^  .  (See 
Figure  8.) 
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Finally,  let  N(r,  s)  be  the  number  of  (r+s )-simplices  this  process 
produces  in  a  partitioning  of  A^_  *  .  Then  clearly 

N(r,  s)  =  N(r-1,  s)  +  N(r,  s-1)  =  N(r,  s)  . 


Since  A„*  A  and  A  x  A.  are  already  simplices, 
Os  r  0 

N(0,  s)  =  N(r,  0)  =  1 


and  since  the  rectangle  A^  x  A^  splits  into  two  triangles  (i.e., 
pyramids  on  two  A^  '  s) 


N(l,  1)  =  N(0,  1)  +  N(l,  0)  =  2  . 


Thus,  the  formula 


N(r,  s) 


r+s  ] 

fr+sl 

.  r  J 

1 

(32) 


holds  when  r+s  equals  1  or  2  .  By  induction  and  the  identity 


u-1 

v 


u-1 
v— 1 


the  formula  follows  for  all  nonnegative  values  of  r  and  s  .  □ 


Finally,  we  note  that  due  to  the  "symmetry"  of  A^_  x  Ag  ,  the  vertex 
U  (referred  to  in  the  proof  above)  can  be  picked  arbitrarily  out  of  the 
set  of  (r+l)(s+l)  vertices  of  Ar  *  Ag  .  Moreover,  we  can  construct 
the  partitioning  of  Ar  *  Ag  in  such  a  way  that  all  simplices  of  the 
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triangulation  have  the  same  two  vertices  in  common.  Indeed,  pick  a 

vertex  U  and  an  (r+s-l)-face  of  the  first  type  and  an  ( r+s-l)-face 

F2  of  the  second  type  such  that  U  ?  F^  U  .  Pick  a  vertex  V  of 

the  (r+s-2)-face  F,  Hf„  ,  which  is  of  type  A  .  x  A  ,  .  Then  we  can 
1  l  r-i  s-l 

use  V  as  the  apex  of  the  pyramids  in  the  inductive  partitioning  of 
F^  and  F2  .  Hence  all  the  simplices  of  the  triangulations  of  F^ 
and  F2  have  V  and  U  as  vertices.  (See  Figure  8.)  Note,  however, 
that  while  U  can  be  picked  arbitrarily,  the  choice  of  V  is  not  com¬ 
pletely  arbitrary,  for  after  V  has  been  chosen  there  are  only  rs 
possible  candidates  for  V  as  U  determines  F^  and  F2  uniquely. 
This  can  be  made  clear  as  follows. 

Let  the  vertices  of  A  and  A  be  denoted  a  ,  a  ,  ...,  a  and 

r  s  0  1  r 

bp,  b^,  ...,  bg  ,  respectively.  Then  the  vertices  of  Af  x  Ag  can  be 

represented  by  the  rectangular  array  {(a^,  b^):  0  i  j;  r  ,  0  _<  j  ^  s 


(V  V 

(ao’ 

V 

(ar  bQ) 

(ar 

V 

(a0*  bs> 


(ar  bs) 


(a,. ,  bQ)  (flj. ,  bx) 


(a^ 


bs> 


(33) 


Without  loss  of  generality,  we  may  assume  that  U  =  (a^,  b^)  .  Then 
clearly  the  vertices  of  F^  are  given  by  the  subarray  given  by  deleting 
the  first  row  of  the  array  (33).  The  vertices  of  F2  are  given  by  de¬ 
leting  the  first  column  of  (33).  Thus  the  vertices  of  F^  O Fj  are 
given  by  the  subarray  obtained  by  deleting  the  first  row  and  column  of 
(33).  This  being  done,  V  may  be  chosen  as  (a^,  bg)  . 

Now,  following  a  suggestion  privately  communicated  to  us  by  M. 
Hazewinkel,  we  may  conveniently  describe  a  procedure  for  triangulating 
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A  <  A^  which  at  the  same  time  gives  a  "geometrical"  explanation  for 
the  formula  (32).  Consider  the  set  of  all  paths  in  the  array  (33)  from 

the  upper  left-hand  corner  to  the  lower  right-hand  corner  which  move 

right  or  down  only.  For  each  such  path,  the  indices  of  the  corresponding 

array  entries  form  a  lexicographically  increasing  sequence  of  length 

r  +  s  +  1  ,  and  the  corresponding  vertices  of  A  *  determine  a 

(r+s ) -simplex.  The  number  of  distinct  paths  of  this  kind  is  readily 

seen  to  be  i  53  j  t  and  the  corresponding  (r+s  )-simplices  all  have  the 

vertices  U  =  (a.,  b.)  and  V  =  (a  ,  b  )  in  common.  They  in  fact 
u  U  r  s 

yield  a  partitioning  of  A^_  *  Ag  as  described  above,  as  we  shall  now 

show.  Any  two  of  the  (r +s )-simplices  corresponding  to  two  distinct 

paths  of  the  above  kind  can  have  at  most  r+s  vertices  in  common  as 

the  corresponding  paths  can  have  at  most  r+s  entries  of  the  array 

in  common.  Hence  all  the  simplices  corresponding  to  the  above  paths 

have  pairwise  disjoint  interiors,  i.e.,  have  at  most  a  common  (r+s-1)- 

face.  To  show  that  they  in  fact  exhaust  A  x  A  we  proceed  as  follows. 

r  s 

Clearly  their  union  is  contained  in  Ar  x  Ag  .  Now  choose  an  arbitrary 
point  (u,  v)  £  4^  *  4g  .  Its  membership  in  one  of  the  (r+s)-simplicies 
can  be  determined  by  applying  the  so-called  Northwest  Corner  Rule  to 
the  transportation  problem  [3,  p. 361]  in  which  the  barycentric  coordi¬ 
nates  of  u  €  A^  and  v  €  Ag  are  the  "supplies"  and  "demands",  respect¬ 
ively.  The  cells  (i.e.,  array  locations)  of  the  basic  variables  so 
chosen  correspond  to  the  vertices  of  a  (r+s)-simplex  (in  the  triangula¬ 
tion  of  A  x  Ag  )  which  contains  (u,  v)  .  Thus,  we  have  a  triangulation 
of  A  x  A  This  is  just  the  partitioning  of  A  *  A  described 
earlier  can  be  seen  as  follows.  All  the  paths  specified  above  which 
leave  U  vertically  (horizontally)  yield  to  a  triangulation  of 
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(F2)  which  is  just  the  inductively-assumed  partitioning  of  (F^) 
described  earlier.  This  provides  a  convenient  scheme  for  determining 
a  triangulation  of  Af  x  Ag  . 


-  BMim 
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